Brightness of a phase-conjugating mirror behind a random medium 
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A random-matrix theory is presented for the reflection of 
light by a disordered medium backed by a phase-conjugating 
mirror. Two regimes are distinguished, depending on the 
relative magnitude of the inverse dwell time of a photon in 
the disordered medium and the frequency shift acquired at 
the mirror. The qualitatively different dependence of the re- 
flectance on the degree of disorder in the two regimes suggests 
a distinctive experimental test for cancellation of phase shifts 
in a random medium. 

PACS numbers: 42 .65.Hw, 42.25.Bs, 42.68.Ay, 78.20.Ci 
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A phase-conjugating mirror has the remarkable ability 
to cancel phase shifts between incident and reflected light 
|^|~^. This cancellation is used in optics to correct for 
wave front distortions Q : A plane wave which has been 
distorted by an inhomogeneous medium is reflected at a 
phase-conjugating mirror; After traversing the medium 
for a second time, the original undistorted wave front is 
recovered. It is as if the reflected wave were the time 
reverse of the incident wave. 

Wave front correction is possible if the distorted wave 
front remains approximately planar, because perfect time 
reversal upon reflection holds only at a single angle of in- 
cidence. For this reason phase-conjugating mirrors have 
been studied mainly in combination with weakly inhomo- 
geneous media. One might think that the diffusive illumi- 
nation resulting from a strongly inhomogeneous medium 
would render the effect of phase conjugation insignifi- 
cant. In this paper we show that, on the contrary, phase 
conjugation drastically modifies the reflectance even if 
propagation through the medium is completely diffusive. 

The phase-conjugating mirror we consider consists of 
a material with a large non-linear susceptibility, which 
is pumped by two counter-propagating beams at fre- 
quency loq. a disordered medium (length L, mean free 
path /) is placed in front of the mirror, and illuminated at 
the other end by a diffusive source at frequency wq + Aw 
(see Fig. |l|) . The reflected light contains both frequencies 
ujQ + Auj and coq — Aw. The reflectances i?_ are de- 
flned as the reflected power at frequency lu± = uq ± Auj 
divided by the incident power. (Note the difference be- 
tween the system we study here and the system studied 
in e.g. Ref. where phase conjugation occurs inside 
the random medium itself.) We distinguish a degenerate 
and a non-degenerate regime, depending on the relative 
magnitude of the frequency shift Auj and the inverse of 



the mean dwell time TdwoU — L"^ /cl of a photon (velocity 
c) in the disordered medium. In the degenerate regime, 
Aw <C l/rdwcU, phase conjugation leads to constructive 
interference of multiply scattered light in the disordered 
medium. In the non-degenerate regime. Aw 3> l/Tdwcii, 
interference effects are insignificant. A distinguishing fea- 
ture of the two regimes is that the reffectance i?_ de- 
creases monotonically as a function of L jl in the degen- 
erate regime, while in the non-degenerate regime it first 
decreases and then increases. The disappearance of the 
refiectance minimum on reducing Aw, provides a distinc- 
tive experimental test for phase-shift cancellation in a 
random medium. 
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FIG. 1. Schematic drawing of the disordered medium 
backed by a phase- conjugating mirror (PCM). Light inci- 
dent at frequency wq -f Aw is reflected at the two frequencies 
Wo ± Aw. 

Our theoretical approach applies techniques of 
random-matrix theory |^] which were originally devel- 
oped for the study of phase conjugation of electrons by a 
superconductor |^ . To reduce the optical problem to the 
scattering of a scalar wave, we choose a two-dimensional 
geometry. The scatterers consist of dielectric rods in the 
^-direction, randomly placed in the x-y-plane. The elec- 
tric field points in the z-direction and varies in the x-y- 
planc only. Two-dimensional scatterers are somewhat 
artificial, but can be realized experimentally We 
believe that our results apply qualitatively to a three- 
dimensional geometry as well, because the randomiza- 
tion of the polarization by the disorder renders the vector 
character of the light insignificant. 

The z-component of the electric field at the frequencies 
w± is given by 



E±(x, y, t) = Y{.cE±{x, y) exp(— iw±f). 



(1) 



The phase-conjugating mirror (at x ~ Q) couples the two 
frequencies via the wave equation [l|J^,p^ 
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The complex dimensionless coupling constant 7 is zero for 
a; < and for x > L,., with Lc the length of the non-linear 
medium forming the phase-conjugating mirror. For < 
X < Lc we put 7 = joe^'^ . The Hclmholtz operator Hq 
at frequency ujq is given by 



(3) 



where ka = uq/c and s(x,y) is the relative dielectric con- 
stant of the medium. We take e = 1 everywhere except in 
the disordered region —L < x < 0, where e = l+5e{x, y). 
The spatial fluctuations Se lead to scattering with mean 
free path I ^ l/fep. The validity of Eq. (||) requires 
Aw/wo ^ 1 and 70 ^ 1. The ratio of these two small 
parameters 



5 = 2Aa;/7oijJo 



(4) 



is arbitrary. 

To define finite-dimensional scattering matrices we em- 
bed the disordered medium in a waveguide (width W)^ 
containing N± — Int(w±W^/c7r) 3> 1 propagating modes 
at frequency ljJ±. A basis of scattering states consists of 
the complex fields 

E±,n = k±]i^ sin (^^) exp(i/c±,„x - iu}±t), (5a) 

^±,« = ^±]i^ sin {-^^ exp{-ik±^nX ~ iLU±t). (5b) 

Here n — 1,2,..., N± is the mode index and the super- 
script > (<) indicates a wave moving to the right (left), 
with frequency uj± and wavenumber fc±^„. The normal- 
ization in Eq. (^ is such that each wave carries the same 
flux. 

With respect to this basis, incoming and outgoing 
waves are decomposed as 



n=l 



n—1 71—1 



(6a) 



(6b) 



The complex coefficients are combined into two vectors 
u = (u+,i,'u+,2, ■ • .,u+^N+,u*_j^,u*_2, ■ ■ -^uI^.nJ^, (7a) 



The reflection matrix r relates u to v. 



V = ru, r = 



r I r — 



(8) 



The submatrices r±± have dimensions N± x N±. The 
reflectances R± are defined by 



R- = N^^ Trr_+rL+, R+ = N^^ Trr+_^r 



(9) 



For Auj <C we may neglect the difference between A''_|_ 
and N- and replace both by iV = lnt{koW/TT). 

We construct r from the reflection matrix rpcM of the 
phase-conjugating mirror and the scattering matrix S of 
the disordered medium. In the absence of disorder, an in- 
coming plane wave in the direction (cos (j), sin </>) is retro- 
reflected at the phase-conjugating mirror in the direction 
(— COS0, — sint/)), with a different frequency and ampli- 
tude. The 2N X 2N reflection matrix rpcM is [|l|,|jlO) 



rpcM 







-lae 




-iip 



(lOa) 



anm = ^nma(</>ri), 4>n arcsin(n7r/fcoM^), (10b) 
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(10c) 



with a = ^7ofcoic- The disordered medium in front of 
the phase-conjugating mirror does not couple uj^ and . 
Its scattering properties at frequency uj are described 
by two N X N transmission matrices t2i(w) and ti2(w) 
(transmission from left to right and from right to left), 
plus two N X N reflection matrices ru (w) and r22 (1^) (re- 
flection from left to left and from right to right). Taken 
together, these four matrices constitute a 2N x 2N scat- 
tering matrix 



\^t2l(w) r22(w) 



(11) 



which is unitary (because of flux conservation) and sym- 
metric (because of time-reversal invariance). (In con- 
trast, rpcM is not flux conserving.) Without loss of gen- 
erality the reflection and transmission matrices of the 
disordered region can be decomposed as 



rii(c^±) - zU±^U^, t2i(cc;±) = V±V^U^, (12a) 
r22{^±) - «V±VpIVL ti2(c^±) = U±0^V^. (12b) 



Here U± and V± are N x N unitary matrices, and 
T± = 1 — is a diagonal matrix with the transmis- 
sion eigenvalues T± „ G [0, 1] on the diagonal. 

Combining Eqs. (joj) (O) we flnd expressions for R± in 
terms of t± and fl = VLaV-(-. The expression for i?_ is 



i?_ = N-^Ti-T-ft (1 - ^/p7 

• (1 - J7^^/prr2*^ 



(13) 



The expression for i?+ is similar (but more lengthy). To 
compute the averages {R±) we have to average over t± 
and V±. We make the isotropy approximation |^ that 
the matrices V± are uniformly distributed over the uni- 
tary group U{N). For TdwciiAo; ^ 1 we may identify 
V+ = V_, while for TdwcuAci; 3> 1 the matrices V-|_ and 
V_ are independent. In each case the average OYevU(N) 
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with A'' ^ 1 can be done using the large A'^-expansion of 
Ref. ||l2|. The remaining average over T± „ can be done 
using the known density p{T) of the transmission eigen- 
values in a disordered medium Q . 

In the non-degenerate regime (rdwcii Acj ^ 1) the result 

is 



{R-) = 



{R+) = 1 - To 



(l-To)2A2' 

Ti{l~n)A^ 



l-(l-To)M2^ 



(14a) 
(14b) 



where To = (1 -I- 2L/Trl)~^ is the transmittance at fre- 
quency loq of the disordered medium in the large- limit 
psf . The quantity A = A^^^Traa^ is the modal aver- 
age of the reflectance of the phase conjugating mirror 



{A J^^^d(l)\ai(t>)\^cos(l) for N 



oo 



Eq. (|T|) 



also be obtained within the framework of radiative trans- 
fer theory, in which interference effects in the disordered 
medium are disregarded 
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FIG. 2. Average reflectances (-R±) as a function of L/l 
for a = 7r/4 and 5 = 0.6, 0.9. The dashed curves are the 
non-degenerate case, given by Eq. ( |l4|). The solid curves are 
the degenerate case, given by Eq. ( [l5[ ) for L/l > 3. Data 
points are results from numerical simulations (open symbols 
for the non-degenerate case, filled symbols for the degenerate 
case) . Error bars are the statistical uncertainty of the average 
over 150 disorder configurations. (When the error bar is not 
shown it is smaller than the size of the marker.) 

In Fig. ^ we have plotted the result (|lj) for (R±) in the 
non-degenerate regime, as a function of L/l for a = 7r/4 
and two choices oi 6 — 0.6, 0.9 (dashed curves). For 
A > 1 (corresponding to (5 = 0.6) the reflectance (i?-) 
has a minimum at L/l = \Tr{A'^ — 1)^^, and both {R+) 
and diverge at L/l = ^7r(y4 - 1)"^ This diver- 

gence is preempted by depletion of the pump beams in 
the phase-conjugating mirror, and signals the breakdown 
of a stationary solution to the scattering problem. For 



A < \ (corresponding to 5 = 0.9) {R-) tends to as L"^ 
for L — > oo, while approaches 1 as 

The situation is entirely different in the degenerate 
regime (rdwcuAcij <C 1). The complete result is a com- 
plicated function of L/l (plotted in Fig. ||, solid curves). 
For L/l ^ oo the result takes the simpler form 



(i?_) = 2roRe^oH. 



1) 



,*2 



arctanhoo, 
-j2^ arctanhoQ, 



"0 "-0 

where the complex number is determined by 

i'Tt/2 







cos(/) a((/)) 
1 - ao a{(j)) 



1 



a, 



(15a) 
(15b) 

(15c) 







When (5 ^ 0, ao ^ 1.284 - 0.0133 i for a = 7r/4. Both 
{R~) and {R+) have a monotonic L-dependence, tending 
to and 1, respectively, as 1/L for L ^ oo. 

To test the analytical predictions of random-matrix 
theory we have carried out numerical simulations. The 
Helmholtz equation 



{\/^ +eujl/c^)S ^ 



(16) 



is discretized on a square lattice (lattice constant d, 
length L, width W). The relative dielectric constant e 
fluctuates from site to site between 1 ± fc. Using the 
method of recursive Green functions [ p^ we compute 
the scattering matrix S{uj) of the disordered medium at 
frequencies uj^ and cij_. The reflection matrix rpcM of 
the phase-conjugating mirror is calculated by discretizing 
Eq. (^). From S(a;±) and rpcM we obtain the reflection 
matrix r of the entire system, and hence the reflectances 
i?± using Eq. (I) 

We took W = bid, Se = 0.5, a — 7r/4, and var- 
ied 6 and L. For the degenerate case we took w+ ~ 
uj- = 1.252 c/d, and for the non-degenerate case 0-"+ = 
1.252 c/d, ~ 1.166 c/d. These parameters correspond 
to 7V+ = 22, Z+ = 15.5 d at frequency uj+. (The mean 
free path is determined from the transmittance of the 
disordered region.) In the non-degenerate case we have 
— 20, I- = 20.1 d. For comparison with the analytic 
theory, where the difference between and and be- 
tween Ijf- and L is neglected, we use the values iV+ and 
l+. Results for the average reflectances (averaged over 
150 realizations of the disorder) are shown in Fig. H, and 
are in good agreement with the analytical predictions. 

A striking feature of the degenerate regime is the ab- 
sence of the minimum in as a function of L/l for 
A > I. A qualitative explanation for the disappear- 
ance of the reflectance minimum goes as follows. To first 
order in L/l, disorder reduces the intensity of light re- 
flected with frequency shift 2Auj, because some light is 
scattered back before it can reach the phase-conjugating 
mirror and undergo a frequency shift. To second order in 
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L/l^ disorder increases the intensity because it traps the 
light near the mirror, where it is amphfied by interaction 
with the pump beams. This explains the initial decrease 
of (i?-) followed by an increase in the non-degenerate 
regime. The decrease persists in the degenerate regime, 
because resonant transmission through the disordered re- 
gion makes trapping inefficient. The resonant transmis- 
sion is the result of constructive interference of multiply 
scattered light, which is made possible by phase-shift can- 
cellation. 

In conclusion, we have studied the interplay of opti- 
cal phase-conjugation and multiple scattering in a ran- 
dom medium. The theoretical predictions of a reflectance 
minimum provides a clear signature for experimental- 
ists in search for effects of phase-shift cancellation in 
strong inhomogeneous media. The random-matrix ap- 
proach presented here is likely to have a broad range of 
applicability, as in the analogous electronic problem . 
One direction for future research is to include a second 
phase-conjugating mirror opposite the first, with a dif- 
ferent phase of the coupling constant. Such a system is 
the optical analogue of a Josephson junction | [To[ |, and it 
would be interesting to see how far the analogy goes. 
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